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The Poisson-Boltzmann (PB) theory is extensively used to gain insight on charged 
colloids and biological systems as well as to elucidate fundamental properties of inter- 
molecular forces. Many works were devoted in the past to study PB related features and 
to confirm them experimentally. In this work we explore the properties of inter-plate 
forces in terms of different boundary conditions. We treat the cases of constant surface 
charge, constant surface potential and mixed boundaries. The interplay between electro- 
static interactions, attractive counter-ions release, and repulsive van 't Hoff contribution 
are discussed separately for each case. Finally, we discuss how the crossover between 
attractive and repulsive interactions for constant surface charge case is influenced by the 
presence of multivalent counter-ions, where it is shown that the range of the attractive 
interaction grows with the valency. 

1. Introduction 

The first applications of Poisson-Boltzmann (PB) theory date back to the beginning of 
the 20th century with the pioneering works of Gouy and Chapman [[H [2], which dealt 
with the counter-ions cloud formed in the vicinity of charged plates. A decade later, 
Debye and Hiickel [ i3| worked out a theory for electrolyte solutions, based on the linear 
PB theory. Since then the PB theory was extensively used for many basic and applied 
purposes, such as the seminal DLVO theory [ HI [5] that explains the stability of colloidal 
suspensions. 

Nowadays, the PB theory is a benchmark tool for a variety of experimental and in- 
dustrial applications [ El |6], in which it is being used to interpret experimental data 
quantitatively. It is a vital ingredient in any attempt to explain intermolecular forces and 
interactions between charged macromolecules, particles and surfaces. 

However, the theory has some drawbacks, which were studied during the last decades. 
Some of the refinements are related to the mean-field limitations of the theory that does 
not account correctly for correlations and fluctuations [HI E]. Other modifications take 
into account additional interactions, such as steric effects [El [10] and hydration interaction 
[ m] . All these refinements are motivated by a number of cases where the PB theory fails 
to explain experimental observations. For example, the theory misses some important 
features when treating multivalent ions and highly charged objects. 

Nevertheless, as the collection of problems that one can treat with the PB theory is 
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very broad in its scope and range, we aim in the present work to further explore several 
PB related cases. We focus on two main aspects. In section [21 the force between two 
charged plates is discussed as function of boundary type. In section [31 we extend the 
discussion to multivalent ions and discuss their effect on the crossover from attractive to 
repulsive pressure. 



2. PB Equation 

Consider a two-plate system occupying the space between two planar surfaces located at 
y=0 and at y=D {y is the coordinate perpendicular to the plates). The system is coupled 
to a bulk reservoir of salt concentration rib, assumed here to be monovalent, z=l. The 
PB equation, describing the electrostatic potential and the ionic distributions in between 
the plates, reads: 

where e is the electron charge, e the solvent dielectric constant, k-Q the Boltzmann's 
constant, and T is the temperature. By defining a dimensionless potential = eip/k-QT, 
the PB equation is written in a dimensionless form: 

: sinh , (2) 



dx^ 

where x=Ky is the dimensionless coordinate, rescaled by the inverse Debye screening 
length defined by K^='&'Ke^n}y/ {ek^T). 

Next, the boundary conditions on the plates can be defined in a general manner. The 
boundary condition at ?/=0 and y=d are given by: 

/o(0(O),0'(O))=O and /<, =0 , (3) 

where d = kD is the dimensionless inter-plate spacing. The functions /o and fd are 
chosen with respect to the surfaces properties. This choice determines the nature of the 
interaction between the plates as is elaborated below. 

Using the solutions of the PB equation one can calculate global thermodynamic quanti- 
ties that are accessible in experiments. In particular, it can be shown [[12] that integrating 
once the PB equation gives the osmotic pressure, 11, in terms of the electrostatic potential 
and its electric field, 0': 

(4) 

Here the osmotic pressure is rescaled by 2kBTnb, yielding the dimensionless osmotic pres- 
sure n. Note that the pressure is homogeneous and has no dependence on the coordinate 
X. It is straightforward to calculate the free energy per unit area F (in a dimensionless 
units) : 

F=- f^dxUix). (5) 
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Although not always accurate, it is instructive to solve the linear PB equation 



(6) 



whose solution is given by (f)=A cosh x + B sinh x . The constants A and B are obtained 
by satisfying the boundary conditions 

foiA,B)=0, (7) 

fd {A cosh d + B sinh d, A sinh d + B cosh d) =0 , (8) 

and are functions of the separation d. The osmotic pressure depends on these two con- 
stants: 

U=-A^ - -B^ . (9) 

In the following we choose three different types of boundaries, and discuss how these 
choices affect the nature of the interaction between the plates. 

2.1. Two constant charge plates 

In many physical setups, the surface charge density can be regarded to a very high 
accuracy as a constant. Therefore, many theoretical works addressed this case in the past 
solving it numerically [[l3l[Tl[T5l[T6l[ni[lHl[19]or analytically [ EQl EH ESI [23l El] . The 
functions /o and fd are given by: 

/o(</.(O),0'(O))=0'(O) + ao (10) 
fd{m,(^'{d))=4>'{d)-ad, (11) 

where ctq and ad are the two rescaled (dimensionless) surface charge densities (Tq — > 
4:7ieao/ (nekBT) and ad — )■ 4:Tcead/{KskBT) , respectively. Substituting them in eqs. [7] 
and [8], we obtain 

A={ad + ctq cosh (i) /sinh (i and B= — ctq • The pressure and the free energy, given by eq. [9] 
and [5], read 



^_ 2aoad cosh d + al + al 

2sinh2c/ ' ^ ' 

2aQad + {al + a^)e-'^ 



2 sinh d 



(13) 



For large separations [d^ 1), the pressure decays exponentially 11 ~ a^ad'^''^ , where the 
sign of the pressure is determined by the relative sign between ctq and ad [EO]. For unlike- 
charges (ctq • ad < the pressure is attractive (11 > 0), while like-charges (ctq ■ o"d > 0) 
repel each other. For small separations there is no dependence on the sign of ctq and 
ad and 11 becomes purely repulsive: n=(cro + 0-^^)^/(2 sinh^ rf) > 0, because of entropic 
reasons (see below). 

Thus, there is a crossover from attractive to repulsive pressure for unlike- charges. There 
is one exception for this finding: when aQ= — ad the pressure is purely attractive. For 
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like-charges, the pressure is monotonically positive and the charges repel each other for 
all the separations. 

The physical interpretation for this behavior is based on the interplay between two 
mechanisms as was previously discussed {e.g., see Ref. [[2l])- When the plates are brought 
closer, the neutralizing ionic clouds increasingly overlap. Since the ionic clouds have 
opposite signs, pairs of negative and positive counter-ions gain additional entropy by 
staying in the bulk instead of neutralizing the system, while the plates, in turn, neutralize 
each other. The closer the plates are, the larger is the ionic cloud overlap, and more 
counter-ions are released to the bulk resulting in an attractive pressure. However, when 
the magnitude of the surface charges is dissimilar {\ao\ ^ a finite amount of counter- 
ions are forced to stay in between the plates to neutralize the excess of charge. The 
concentration of these neutralizing counter-ions grow when the separation is reduced, 
leading to a repulsive pressure (due to the van 't Hoff pressure of the counter-ions). For 
small separations (rf — )■ 0), this repulsion diverges as 11 ~ l/c?. 

For the pure antisymmetric case, where oq= — Od , only the counter-ion release mech- 
anism is observed since no counter-ions are needed to neutralize the plates in small sep- 
arations (there is no charge excess and the two plates completely neutralize each other). 
This is the origin of the pure attractive pressure in this case. On the other hand when 
o'o ■ > 0, only the repulsive pressure, caused by counter- ion neutralization, takes place 
since the clouds have similar signs, and there is no counter-ions release. We hence con- 
clude that the crossover from attractive to repulsive pressure occurs only for unlike-charges 
(ctq • (Trf < 0) with asymmetric magnitudes, \oq\ ^ \ad\. 

In Figure la three curves of pressure vs. distance with constant surface charge bound- 
aries are shown. In addition to the simple analytical results obtained for the linearized PB 
case, we solve numerically the non-linear PB equation. Both are presented in Figure 1; 
the linear case as dashed lines while the non-linear solid lines. The comparison 

between the solid and dashed curves shows that the solution of the linear equation agrees 
qualitatively with the exact solution of the non-linear equation. The numerical deviations 
between the two are small for the chosen values of ctq and a^. However, for larger surface 
charges the deviations grow, and the linear solution is not sufficiently accurate anymore. 



2.2. Two constant surface potentials 

The functions /o and fd are given here by: 

/o(0(O),0'(O))=0(O)-0o 
/,(0(rf),0'(rf))=0(rf)-0rf. 



(14) 
(15) 



where 0o and 0^ are dimensionless potentials. Obtaining the constants A and B from these 
two relations, v4=0o and B={(f)d — 0o cosh ci) /sinh ci , leads to the dimensionless osmotic 
pressure H and free energy F: 



20o0d cosh d 



2 sinh^ d 



20, 



09d 



^1 + 



2 sinh d 



(16) 
(17) 




Figure 1. Osmotic pressure 11 as function of inter-plate separation, d. The solution 
to the linear (dashed) and non-linear (solid) PB equation are presented for constant 
surface charge case in (a) and constant surface potential in (b). Three profiles are 
shown for each case: I, II and II indicate the repulsive, attractive and a crossover pro- 
file, respectively, (a) The boundary conditions for profiles I, II and II are: ao=ad=S, 
(To=— cr(/=3 and (7o=—3ad/2=3, respectively, (b) The boundary conditions for I, II 
and II are: ^o=?^d=3, ^o=— ?^d=3 and ^o=— 3^(^/2=3, respectively. 



These expressions are almost identical to the constant surface charge case (except for 
the sign), and indeed, for large separations the relative sign of 0o and (f)d determines the 
pressure: IT ~ ^o^^e"'^ , similarly to the large separation behavior of the constant charge 
case. However, for small separations, the pressure is different: 11^ — (0o "~ 0d)^/2sinh^ d, 
yielding a pure attractive pressure, except for 11=0 when 4>q=4>ci- Unlike the constant 
charge case, in the constant potential case the counter-ion concentration remains constant 
near each plate, but the effective surface charge diverges when the plates are brought closer 
together. This results in a diverging electrostatic attraction that governs the pressure for 
small d. 

Another remark about the constant potential case is that the crossover from repulsive 
to attractive pressure prevails only for like-potentials (0o " 0(i>O). For unlike-potentials 
(00 ■ 0<i<O) the pressure is purely attractive, while for the special symmetric case, (f)o—(f)d , 
the pressure is purely repulsive. 

In Figure lb we show three curves of pressure vs. distance for three different surface 
potential values (attractive, repulsive and crossover). Note that the attractive profile 
diverges, while the repulsive one saturates and is different than the constant surface 
charge case, where the repulsive one diverges and the attractive one saturates. This is 
due to the different constraints that the boundaries apply on the ionic solution. Also 
note that the dominant deviations of the exact non-linear PB sohition from the linearized 
one are large for small inter-plate separations, due to the large overlap between the two 
induced counter-ion clouds. 
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2.3. The mixed case 

The third case is a mixed boundary condition, where the plate at y=0 has a constant 
potential while the one at y=D has a constant charge. The functions /o and are 
given by: /o (0(0), 0'(O)) =0(0) - 0o , and fa {(f){d) , (j)' (d)) =(j)'{d) - • Substituting the 
solution for (p in these two relations, the constants A and B are obtained: A=(f)Q and 
B={afi — 00 sinh d) /cosh (i . The pressure 11 and the free energy F are given by: 

^_ 20oarfsinhrf + 0g-org ^^^^ 
2 cosh^ d ' 

^_20ofTd + (0g-a2)e_2^ ^^g^ 



2 cosh d 

The large separation the pressure behaves as n~0oO"de~'^ , similarly to the previous two 
cases where the nature of the pressure is determined by the relative sign (0o ■ o^d)- The 
similarity between the three cases for large d is due to the decoupling between the plates, 
effectively leading to two isolated planes with no mutual effect on their surface potential 
or charge. On the other hand, examining the small separation limit, we obtain 
n~(0Q — (j^)/(2cosh^ d) where the sign of the pressure depends on the magnitudes of 
00 and (Trf, unlike the former cases. When |(Trf|>|0o| the pressure is attractive, while in 
the opposite case, |(jrf|<|0o|, the plates repel each other. The competition between the 
counter-ions release, the charge neutralization, and constant ion density at the constant 
potential surface leads to this versatile behavior. Moreover, in contrast to the former 
cases, there is no divergence here at small separations, and the pressure saturates as 
(i — )■ 0. Note the special case where |c'"rf| = |0o| that leads to vanishing of the pressure as 
d^Q. 

In this section we demonstrate how versatile can be the electrostatic pressure between 
two plates in the framework of the linear PB theory. We note that the three cases treated 
here have the simplest boundary conditions one can consider. In physical realizations, 
one should often take also into account the charge regulation mechanisms [ 125j , where 
the surface properties are determined by additional surface parameters such as the ionic 
dissociation degree. In these cases the functions /o and fd have a more complex form, that 
describes the relation between and 0' at x=0 or x=d due to surface activity. Although 
our results are derived for the linearized PB case, we expect similar conclusions to hold 
for the general non-linear PB case as well. 

3. Attraction to Repulsion Crossover with Multivalent Salts 

We consider a two-plate system with unlike- charges, o"o<0 and crrf>0 in the presence of 
ionic solution with multivalent salt. First we treat the case of a symmetric salt {z:z) and 
then consider mixtures of two salts: a symmetric monovalent salt (1:1) combined with an 
asymmetric multivalent salt (1:^;). 

3.1. Multivalent symmetric salt {z:z) 

In this case the dimensionless PB equation of eq. [2] is generalized and reads: 

-4=2sinh(z0). (20) 
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where z is the ionic valency. The boundary conditions at x=0 and d are given by 

/o(</.(O),0'(O))=0'(O) + aoA and /, =0'(rf) - a,/^ . (21) 

In this form the equation is almost identical to the monovalent salt case (eq. HJ. Therefore, 
the attraction to repulsion condition of the monovalent case (as was discussed in Ref. [ 
EH) can be used: 



Yi=l <^ e-v^*^ < 7o/7d < e^"^ and (Tq ■ < ,^2) 
' > otherwise 



where 7o,d=\/ {2y/z/ao^df + 1 - \2y/z/ao^d\- 

3.2. Mixtures of electrolytes 

Next we consider a binary salt mixture: the first is an asymmetric salt (l:^;) for which 
the bulk concentration of the multivalent cations and monovalent anions are nl and zfif , 
respectively. The second is a monovalent salt (1:1) for which the bulk concentration . 
The dimensionless PB equation now reads: 

za 



sinh0- — (e-"^-e^) . (23) 



where a=nl / Ub is the ratio between the bulk densities of the multivalent and monovalent 
cations. This equation depends explicitly on the parameters a and z, while all the other 
physical parameters are taken into account by the dimensionless boundary conditions 
(eqs. [To] and [TT!) . The first integration of this equation gives the dimensionless pressure 
in between the plates: 

ni„= -1<P" + cosh + |e-^* + (24) 
The pressure in the outer region is calculated by taking (p=(p^=0: 

nout=i + 2 + f ^^^^ 

In order to calculate the attraction/repulaion crossover line we examine the point in the 
parameter space where the net pressure on the plates 11=11^ — Ilout vanishes, yielding the 
following relation between and 0': 



2(cosh - 1) + a(e-^'^ - 1) + za{e^ - 1) (26) 

By solving this first-order differential equation, the boundary line between attraction and 
repulsion regions in the {d, |o"rf/cro|) plane can be obtained. 

In the limit of low salt densities, one can assume that the concentration of the mul- 
tivalent counterions is much larger than that of the monovalent ones due to entropy 
considerations. Hence, we can solve the counter-ions only equation, and derive a similar 
condition for low salt limit as was derived in Ref. 

1 1 



< A ' (27) 
a* 
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Figure 2. Crossover between attraction and repulsion for binary salt mixtures, (a) 
Different values of the multivalent cation concentration (see box). The monovalent 
salt concentration is n5=10~^ M and the negative surface charge density is (To= — 
e/lOOA^^. The valency of the cations is z= + 3. (b) Varying the multivalent cation 
valency, z. The negative surface charge density ao is the same as in (a). The overall 
ion bulk concentration, nt=2nb + {1 + z)ani,=4nb , is kept constant for all the curves, 
yielding a constant reference pressure in the outer region. 



where a*=ekBT/ {27ieD) is a characteristic charge surface density that depends on D. 
This condition has no dependence on the salt concentration, so that the multivalent and 
monovalent cases have the same behavior in this limit. 

In Figure 2a we show the effect of adding multivalent salt to the solution. It is evident 
from the plot that the attraction region grows with adding multivalent salt. In Figure 2b 
the effect of increasing the valency of the multivalent cations is examined. In the plot it 
can be seen that the attraction region is increased for higher valencies. It is interesting to 
note that the crossover distance d scales inversely with the valency for constant asymmetry 
ratio |(7(i/a"o|. 

4. Summary 

In this work we explored a variety of PB theory results. In particular, wc focused 
on several types of charge boundaries that affect the nature of the pressure between 
two plates, as well as how multivalent ions influence the crossover between attractive to 
repulsive inter-plate pressure. 

It is evident from this work that solutions of the linear PB equation can serve as a 
qualitative tool for understanding the interaction. However, when interpreting experi- 
mental data, one should solve the non-linear PB theory in order to have better numerical 
accuracy, especially for small inter-plate separation, d. 

Finally, we would hke to stress that the surfaces cannot always be regarded as having 
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either a constant surface charge or potential. In quite a number of experimental setups 
the surface properties depend explicitly on the inter-plate separation. For example, the 
degree of dissociation of surface charge groups can vary with d. In order to take this effect 
into account one should use boundary conditions that are determined self-consistently as 
function of separation [ |25] . 
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